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' Abstract. We compute the eigenfunctions and eigenvalues of the periodic integrable spin s 

XXX model using the coordinate Bethe ansatz. To do so, we compute explicitly the Hamil- 
tonian of the model. These results generalize what has been obtained for spin \ and spin 1 
, chains. 
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CN ■ 1 Introduction 

> 



The resolution of Heisenberg spin chain pQ was initiated in H. Bethe's seminal paper [2] where 
he used a method called now coordinate Bethe ansatz. Since this work, several new methods 
appeared: algebraic Bethe ansatz OH], functional Bethe ansatz (or separation of variables) [5] 
ON ■ or analytical Bethe ansatz [6]. These more elaborated techniques allowed one to go further: new 

solvable models have been discovered and new results have been computed such as correlation 
functions. As a consequence, the coordinate Bethe ansatz was neglected. However, this method 
is the simplest one and gives a very efficient way to construct explicitly eigenfunctions, but it 
is believed that it works only for simple models. In this note, we show that actually it can be 
applied also to more complicated models as the spin s XXX model. 

This paper is organized as follows. In Section [21 we compute the Hamiltonian of the spin s 
chain we want to solve. To our knowledge, the explicit form of the entries of the Hamiltonian are 
written for the first time. We also compte the su(2) symmetry algebra and the pseudo-vacuum, 
a particular (reference) eigenstate. In Section [3J we present the coordinate Bethe ansatz and 
get the Bethe equations obtained previously by the algebraic or analytical Bethe ansatz. We 
conclude, in Section [H on the advantages of this method and on open problems. 



2 Integrable periodic spin s chain 
2.1 Hamiltonian of the spin s chain 

The Hamiltonian of the periodic integrable spin s chain has been computed in [7j thanks to 
a fusion procedure. This Hamiltonian has been expressed as a polynomial of the invariant 
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of su(2) (see for example also [8] for a review). For our purpose, we need to give an explicit 
expression of the Hamiltonian entries. Namely, the Hamiltonian is the following matrix acting 
on (C 2s+1 )® L 

L 

with the periodic condition L + 1 = 1 and the subscript + 1) stands for the two spaces 
where the (2s + l) 2 x (2s + l) 2 -matrix h acts non-trivially. We choose to enumerate the basis 
of C 2s+1 as follows: \s), \s — 1), . . . , | — s), where \m) = \s,m) denotes^ the spin s state with 
s 2 -component equals to m. 

The non- vanishing entries of the matrix h may be parameterized by three integer parameters 

mi, 777-2, n 

2s min(m2,2s— mi) 

h= ^2 Pmi,m2 I s — m l — n )( s — m ll ® I s — m 2 + n){s — mi]. 

mi,m2=0 n=— min(mi,2s— ma) 

From the results found in [7], one may prove that /3 m m2 can be factorised as: 
(_l)n-l 



/3 n 

rmi 



77 



\ 



Ml + 77 \ / M 2 

Mi U 



2s - Mi \ / 2s - M 2 + n 

77 / \ 77 



for 77 > 0, (2.2) 



where we have introduced the notation 

Mi = min(?77i, 2s — 7772) and M2 = min(?772, 2s — mi). (2-3) 

The remaining /3's are given by the relations: 

Mi-1 -. M 2 -i . 
0° = - V — V — 



2s -4 ^ 2s -T 

e=o £=0 



o-ri on 



m 2 ■ 



2.2 su(2) symmetry 

At each site, we have a spin s representation, and the expression of the su(2) generators in this 
representation reads: 

s 

s~ = ^2 \/ (s + n)(s - 77 + l)\n - l)(nj, 

n= — s+1 

8-1 

s+ = ^ v / ( s - re )( s + n + 1 )h + 1 )( n l ) 

n=— s 
s 

s z = n\n){n\. 

n = — s 

x Let us stress that since all the sites have spin s, we do not mention it, and write only the s~ value of the 
states. 
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They obey 



2s z , [s z ,s ± ] = ±s ± , C2 = ( s -) 2 + I( s + s - +s - s +) = s ( s + i) 



We will note s", a = z,± and j = 1, . . . , L, the generators acting on site j. Let us stress that 
these local operators do not commute with the Hamiltonian % given in (|2.ip . However, there is 
a global su(2) symmetry. The generators of this su(2) symmetry take the form 

5* = 5>J and S± = J2sf. 

3=1 3=1 

They obey the su(2) commutation relations [iS 2 ^ 1 * 1 ] = itS 1 * 1 , [<S + ,5~] = 2S Z . Remark that the 
Casimir operator C 2 = (<S Z ) 2 + t;(S + S~ + S~S + ), although central, is not proportional to the 
identity, since we are considering the tensor product of L spin s representations, a reducible 
representation. 

It is a simple calculation to show that [s° +s" +1 , fyjj+i] =0, a = z, ±. It amounts to check 
the following recursion relations on the coefficients m2 : 



,m,2 



V(mi + l)(2s - mi)/3" 1+ljm2 = y^s + n - m 2 + l)(m 2 - n)Ct 

- V / (m 2 + l)(2 S -m 2 )/3^+ 1 m2+1 + V(2a - n - mi)(n + mi + 1)^ 



m 2 ' 

mi ,rri2 



\/mi(2s -m 1 + ^)Pmi-i,m,2 = V ( 2s + n ~ m 2 )(m2 -n + l)/3, 

- Vm 2 (2s-m 2 + l)/3^7 i 1 m2 _ 1 + ^(2* - n - mi + l)(n + m^^^. 

Hence, S a commutes with %. 

Due to this su{2) symmetry, the wave functions can be characterized by their energy, their 
spin and their S z component. In other words, one can diagonalize the Hamiltonian in a sector 
where the operators S z has a fixed value S z = Ls — m. This is done in the next section. 



2.3 Pseudo- vacuum and pseudo-excitations 

We wish to present the construction of the Hamiltonian eigenfunction in the framework of 
coordinate Bethe ansatz for spin s. The spin s = \ case is the Heisenberg chain, solved in [2]. 
It gave the name to the method. For the case s = 1, the method has been generalized in [9]. 

The first step of the coordinate Bethe ansatz [2] consists in finding a particular eigenvector, 
called the pseudo-vacuum, for the Hamiltonian. It is usually chosen as the vector with the 
highest spin. In the present case, it is the unique vector in the sector S z = Ls: 

\0) = \s) <g) - ® \s). 

Using the explicit forms for the /3's, we get h\i \s) (g> \s) = 0. Thus, \0) is a "^-eigenvector with 
vanishing eigenvalue. 

The second step consists in adding pseudo-excitations. These pseudo-excitations are not 
physical excitations (hence the name pseudo-excitation). They are obtained by acting with a 
creation operator ej , conjugated to sj, on the pseudo- vacuum \0). Let us remark that this 
operator in this finite representation does not satisfy (e~) 2 = but rather (e~) 2s+1 = 0. This 
explains the supplementary difficulties to deal with s > ^- Indeed, in the case s = |, no more 
than one pseudo-excitation can be at the same site: we have strict exclusion. In the general case 
of spin s, we have a weaker exclusion. More precisely, we can have up to 2s pseudo-excitations 
at the same site. This behavior appears already for s = 1. 



4 



N. Crampe, E. Ragoucy and L. Alonzi 



3 Coordinate Bethe ansatz for general spin s 

We define a state in the sector S z = Ls — m, for 1 < x\ < x 2 < • • • < x m < L 

\xi,x 2 ,...,x m ) = e~e~ • • -e- m |0>, (3.1) 
where e~ is conjugated to s~: 



E s + n . , - -1 ■ , (2s)! , v / i 
\ — r n-l)(n =55 5 with 5 = > n)(n . 

V s - n + 1 [s — n)\ 

n=— s+1 n=—s v ' 

This choice for e~ is for later convenience. The set of such non- vanishing vectors (i.e. Xj+ 2s > Xj 
for 1 < j < L — 2s) provides a basis for this sector. As already noticed, contrarily to the usual 
case (s = g), several (up to 2s) pseudo-excitations at the same site are allowed, that is to say, 
some Xj's can be equal. The restriction that no more than 2s particles are on the same site 
is implemented by the fact that (e~) 2s = 0. Using the explicit form of e~ , we can rewrite the 
excited states as follows: 

\xi,x 2 , • • • , x m ) = a mi ■ ■ ■ a mk \s) <g> • • • ® \s) (g)|s - mi) ® |s) ® • • • <g> |s) <g>|s - m 2 ) <g> ■ ■ ■ , 

V » ' ^ v ' 

XI— 1 

where m,- is the number of times a;,- appears and 



a. 



2s 
m 



where ( ^ J = z ^ is the binomial coefficient. Let us remark that, if rrij > 2s the 

vector \s — rrij) has no meaning but the normalization a mj vanishes. 

Any eigenvector in the sector S z = Ls — m is a linear combination of the vectors (|3.ip . Then, 
let us introduce the vector 



, X 2 , • • • , Xm) \ x\ , , • • • , Xm) , 

Xl<X2<-<2:m 

where o(xi, X2, ■ ■ ■ , x m ) are complex- valued functions to be determined. As in the case of s = sj, 
we assume a plane wave decomposition for these functions (Bethe ansatz) 

Oi{x\ , • • • , x m ) — 

Ap(k) exp {i(kp\x\ -\ \-k Pm x m )}, 

-Pe6 m 

where & m is the permutation group of m elements and Ap(h) are functions on the symmetric 
group algebra depending on some parameters k which will be specified belovJE Using the fact 
that the states \x±,x 2 , • • • , x m ) form a basis, we can project the eigenvalue equation 

m> m = E^ m (3.2) 

on these different basis vectors to determine the Ap(k) parameters. 

Since H is a sum of operators acting on two neighbouring sites only, one has to single out 
the cases where the x's obey the following constraints: 



2 In the following, to lighten the presentation, the k-dependence will not be written explicitly. 



Coordinate Bethe Ansatz for Spin s XXX Model 



5 



• all the Xj's are far away one from each other and are not on the boundary sites 1 and L 
(this case will be called generic), 

• Xj + 1 = Xj+\ for some j, 

• xj = Xj + \ for some j, 

• Xj = Xj + i = ■ ■ ■ = Xj +mi and Xj + 1 = Xj +mi+ i = ■ ■ ■ = Xj +mi +i+m 2 f° r some positive 
integers mi and m2, 

• x\ = 1, or x m = L. 

As the eigenvalue problem is a linear problem, it is enough to treat the cases where at most one 
of the particular cases appears: more complicated cases just appear as superposition of these 
'simple' cases. 

Projection on \x\, X2, • • • , x m ) with Xj + 1 < Xj+i, Vj, a?i > 1 and x m < L. As 

usual, we start by projecting (|3.2p on a generic vector \x±,X2, ■ ■ ■ ,x m ). This leads to 



Pee m V 7=1 / 



;X m )) = 



which must be true for any choice of generic x's. Therefore, we get for the energy (using the 
explicit forms of the /3's given in Section [2|) 

1 m 

E = V (2-e^ -e~ ik A. 

2s ^ v ' 



After the change of variable 

(3.3) 



jkj = + is 
Xj — is ' 



the energy becomes 
2s 



Xj + s 2 ' 
J= i j 

This form for the energy is the one obtained by algebraic Bethe ansatz [10] . 

Projection on \x\, X2, • • • , x m ) with Xj + 1 = a^j+i for some j. Let us consider now 
the projection of (|3.2j) when two pseudo-excitations are nearest neighbours. Using the form of 
the energy previously found, we get 

£ A P (p ltl - (3% - $} A + GS^aa - fr^e*"' + (/3 X >2 - AliK^^) 
Pe6 m 

x e j(...fcpjX j +fcp (j+1) (H-x- j )... ) _ Q_ 

This equation is trivially satisfied for s > \ since using explicit values we find f>\ 1 — — /3q x = 0, 
/^ri a 2 ~~ P10 = an d Pl,i a 2 — Pq \ = 0. For the case s = |, we find a constraint between 
and ^4pt (where Tj is the transposition of j and j + 1). Explicitly, it is given by (13.5P with 

S = ^ . . . 

Projection on \x\, X2, • . • , a?m) with = Xj+i for some j. For s > k, we must also 

consider the case when several particles are on the same site. Defining Sj, the shift operator 

adding 1 to the j th variable, we get the following relation, when two particles are on the same 

site 

(1 + S- 1 S^ 1 )(S i S i+1 + (2s - l)Si - (2s + l)S i+1 + l)a(..., Xi ,x h ...) = 



2s(l-2s) 
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=> (S i S i+1 + (2s-l)S i -(2s + l)S l+1 + l)a(...,x i ,x i ,...) = 0. (3.4) 
Using the plane waves decomposition, we get the following constraint 

A PTj =a(e ikp i,e ikp (i+v)A P , (3.5) 
where Tj is the transposition of j and j + 1, and we have introduced the scattering matrix 
uv + (2s-l)u-(2s + l)v + l 

<J(U,V) = ; : ; r . 3.0 

v ' ; uv + (2s - l)v - (2s + l)u + 1 v ; 

As in the case s = |, relation (|3.5p allows us to express all the Ap's in terms of only one, 
for instance (where Id is the identity of & m ). More precisely, one expresses P E & m as 
a product of T, and then uses (|3.5p recursively to express in terms of A^. At this point, 
one must take into account that the expression of P in terms of T is not unique, because of the 
relations 

Tf = Id, [Tj,Ti\ = (|j - i\ > 1) and TT m T = T i+1 TT +1 . 

Therefore, for the construction to be consistent, the function a has to satisfy the relations 

a(u, v)a(v, u) = 1, [c(«, v), cr(w, z)\ = 0, 
a(u 1 ,u 2 )a(ui,u 3 )a(u2,u 3 ) = a(u 2 ,U3)a(u 1 ,u 3 )a(ui,u 2 ). 



By direct computation, we can show that (I3.6P indeed satisfies these relations. We can solve 
the recursive defining relations for Ap and we find, with a particular choice of normalisation 
for A^, the following explicit form, for any P E <& m 

A P= T\(^- {eit "; 1){e T~ 1] 

1 1 \ 2s e lkp i — e lkpk 
j<k v 

This scattering matrix becomes after the change of variables u = and v = 
A — fx — i 



cr(A,/x) 



A — fx + i 



Let us remark that after this change of variables the scattering matrix er(A, Li) does not depend 
on the value of spin and is similar to the one obtained for s = |. 

Projection on |a?i, . . . , xi, . . . , Xj, ccj + 1, . . . , Xi + 1, . . . , x m ). One can compute 

mi m.2 

-fmi ,ni2 ("Si? • • • j S'i+mi — 1 > <Si+mi , • • • , 5i+mi+m2 — 1 ) 

+ fm 2 ,mi (5 , j +mi+m2 _ 1 , . . . , S i+mi _ v . . . , S { / 

x a(. . . , Xi, . . . , Xi, Xi + 1, . . . , Xi + 1, . . . ) = 0, (3.7) 

where 

111-2 



P'tni,m2^y^ ^) — ^ ^n*m2- n^mi fto,m2 ^™2 — ' ' ' ^r«2 
n=l 
Ml 



n=l 
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^ 2 + Pl,0f jj ~ «miam 2 2_, ^ 3 j 



and we have used the notation (12.31) and 



y = (2/1,2/2, • • • ,2/mi) witn Vk = Si+k-l, l<k<mi, 
z = (zi,z 2 , ■ ■ ■ , z m2 ) with z fc = 5 i+mi+fc _i, 1 < k < m 2 . 



Relation (|3.7p is implied by f|3.4|) . The sketch of the proof goes as follows. 

First, one can check directly thatH Po,2(Sj, Sj+\) = 0, V j, since Po^{Sj, Sj + \) corresponds to 
relation (|3.4p . The same is true for P2,o(S~ +1 , Sy 1 ) (after multiplication by Sj'^Sj' 



Next, we rewrite (|3.7p as 

-Pmi,m2 (y! z ) "I - -Pm2,mi ( z j y) = 0, (3. 

where we have defined 

y = (2/mi,---,2/2 _1 ,2/r 1 ), i- e - Vk = S^ mi _ k , l<k< mi , 

~\ zi 1 ), i-e. z k = S7^ mi+m2 _ k , l<k< m 2 . 



z 



These variables are such that if Po^{zj, -Zj+i) = Vj, then we have also Po,2(^j,^j+i) = 0. 
Hence, a property valid for P mi/m2 (y;z) will be also valid for P m2jrni (z; y). 
We first focus on 

m m 

Po,m{z) = 2J ana m -n/3o,m z m-n+l ■ ■ • Z m - a m f3^ ^ Zj + a m (/?o,m ~ m /3l,o) ■ 
n=l j=l 

If we suppose that we have variables Zj such that Po2(%', ^7+1) = 0, Vj, then from expres- 
sion (|2.2p . and after some calculation, one can show that 



in 



ziz 2 • • • z m = 1 - m + 2^ Xj zj, 
3=1 



fc=i 



Thus, the polynomial Po,m( z ) becomes a linear function of the ZjS. Looking at the coefficient 
of Zj and at the constant term, one checks that they identically vanish, so that Po,m( z ) = 0. 
Looking at the general polynomial P mi)m2 (y,z) and using the relation 

Q — 7t QTl 

one can rewrite it as 

Pm\ ,m 2 

(y, z) = a mi P ,m 2 ( z ) + OLM2 PoMliVmi-Mr+l, ■•• ,2/mi) + «Mi «M 2 Pm 1 ,m 2 (y), 



mi-Afi / fl o 

E 



R mi ,m2 (y) — — ^2 ( — 1 — 9 ~ + PlflVj ) + 9 (/^0,M 2 ~~ ^Afi,0 + $),mi — /3o,m 2 ) ■ 



'Multiplication on the right by a(. . . , Xj, . . . , Xi, x, + 1, . . . , Xi + 1, . . . ) will be understood during the proof. 
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Thus, to prove relation (|3.8j) . it is enough to show that 

Rm x ,m 2 {y) + flffl2,mi(z) = 0. (3.9) 
Using the expression of Mi and M2, see (|2.3p . it is easy to see that 

mi — M\ = n\2 — M2 = mi2 > 0. 

Two cases have to be distinguished: mi2 = or mi2 > 0. In the first case, relation (|3.9p is 
trivially satisfied. In the second case, equation (j3.9j) can be rewritten as 

ttll2 

which is obeyed if 

Si+j-i + S i+2s+ j_i = 2. 
To prove this last relation, we use recursively (|3.4[) to show 



S k S k+e + i+(?i-ijs k -(j + 1) S k+i = V j. 

Taking k = i + j — 1, I = 2s and using mi + m2 = 2s + mi2 gives the result. 

Hence, relation (j3.8[) is satisfied if the variables are such that Po,2{zj, z j+i) = Po,2(yj, Uj+i) = 
0, Vj. This ends the proof. 

This step concludes the bulk part of the problem, the other possible equations being fulfilled 
by linearity. It remains to take into account the periodic boundary condition. It is done through 
the following projection. 

Projection on |1, • • • , x m ). As usual, this leads to a constraint on the parameters kj. 
It is not surprising since these parameters can be interpreted as momenta: we are quantizing 
them since we are on a line with (periodic) boundary conditions. Namely, this leads to 

Ap ( ex.p(i(k P2 x 2 H h k Pm x m )) 



P&Gr 



exp(i(kpix 2 + k P2 x 3 H h k P{m _i)X m + k Pm L))) = 0. 



Now, we first perform the change of variable in the summation P — > PT\ ■ ■ ■ T m _i in the second 
term of the previous relation. Then, using recursively relation (I3.5P and projecting on inde- 
pendent exponential functions, we get the quantization of the momenta via the so-called Bethe 
equations 

e iLk 3 = Y[a(e ike ,e ik i) for j = l,2,...,m. (3.10) 

Since these equations express the periodicity of the chain, they are equivalent to the ones ob- 
tained through projection on \xi, . . . ,x m -i,L) (as it can be checked explicitly). Thus, we do 
not have any new independent equations through projections, and the eigenvalue problem has 
been solved (up to the resolution of the Bethe equations). 

Note that using the change of variables (|3.3p and the expression (|3.6p for the scattering 
matrix, equations (I3.10P can be rewritten as 

Xj + is yj \( — Xj — i 



n 



One recognizes the usual Bethe equations of the spin s chain [101 [TT] . 
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Action of su(2) generators 

Since the su(2) generators commute with the Hamiltonian, from any eigenfunction \I/ m , one 
can construct (possibly) new eigenfunctions by application of S a , a = z, ± on As already 
mentioned, it is a straightforward calculation to check that 

S z * m = (Ls-m)^ m . 

Moreover, it is part of the ansatz to suppose that the eigenvector ^ m is a highest weight vector 
of the su(2) symmetry algebra, 



Let us stress that for *$> m to be an eigenvector, one has to assume that the rapidities Xj have to 
obey the Bethe equations. In the same way, \F m is a highest weight vector only when the Bethe 
equations are fulfilled. In the context of coordinate Bethe ansatz, there exists no general proof 
(for generic spin s) of it (at least to our knowledge). Note however that for spin |, the proof 
was given in [12]. Nevertheless, one can check the highest weight property on different cases, 
and we illustrate it below by the calculation of S +i &2 an d tS 4 "^. We also show on the 

last example where the proof used by Gaudin does not work anymore for spin s > ^ . 

The fy m vectors should be also related to the ones obtained through algebraic Bethe ansatz 
(ABA). Such a correspondence, for the case of spin ^, has been done in [13] using an iteration 
trick based on the comultiplication [HI [15]. Let us note that in [13] they used the relation 
(T12) 2 = which is not true anymore for s > ^. Their proof must be generalized to apply in 
our case. Let us also notice the other method using the Drinfel'd twist [16]. Moreover, since it 
is known that the ABA construction leads to su{2) highest weight vectors, and assuming the 
same property for the coordinate Bethe approach, it is clear that the two methods should lead 
to the same vectors, up to a normalisation. 

For instance, considering its ABA "counterpart" takes the form 

$1 = XX^AO • • • Tt'H^iX^T^iX,) ■ ■ ■ 2^010), 

x=l 

where T^\X) is the representation of the monodromy matrix at site j: 

T^(X) = -±-( X+ l 5 i x ^. z ) 
X — is\ isj X-iSj ) 

and Ai is the Bethe parameter. This leads to 

i ( X\ + is\ x _. . 

Xi + is ^\X X - is J 

x=l 

that has to be compared with 

L L 

y l = Y^e iklX \x) = Y j e iklX e~ \0). 

X = l 2 = 1 

Using the change of variable (|3.3|) . it is clear that, apart from a normalisation factor, the two 
vectors are equal. 

Calculation of «S~'~\E , i and S^~^2- A direct calculation leads to 

= 2s— ^— (1 -y L )\0) with y = e ikl , 
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which is identically zero using the Bethe equation y = 1. Hence, ^\ is indeed a highest weight 
vector for the su(2) symmetry. 

In the same way, one can compute 



pg6 2 



J^(2a - l)e l{kl+k2)x \x) + Yl 2se l( - kpiXl+kp * X2) {\x 1 ) + \x 2 )) \ . 

x=l l<xi<X2<L 

Using the relation 

A Tl (ki,k 2 ) = a(ki,k 2 )Ai d (k 1 ,k 2 ) 
and the normalisation Ai^(ki,k 2 ) = 1, one gets 

5 + ^ 2 = ^2<(2s- l)( yi y 2 ) x (l + a( yi ,y 2 )) + 2s 

x=l I 



^2 Vl x 

1 ^1 

1-2/2 



+ ^(2/1,2/2) — : y 2 + -j y 2 + cr(yi,y 2 )- y\ 



1-yi I-2/1 

where ?/j = e* fcj , j = 1, 2. 

Now, from the Bethe equations 

2/f = 0"(2/i,2/2) and 2/2 = cr (2/2, 2/i), 
one simplifies it as 

5 + ^ 2 = ^(2/12/2)" {2s - 1 + 2a f-^- 

^ I Vl-?/2 



1-2/2 



1 



+ cr(yi,y 2 ) 



2s - 1 + 2s 



2/2 1-2/1 
2/1 1 



l-2/i 1-2/2 



Finally, the form of the scattering matrix a ensures that the quantity within brackets {• • • } 
vanishes. 

Calculation of <S"'"\E r 3. Performing the same kind of calculation on ^3, we get 



Pee 3 



+ 2sa 23 ai 3 



x=l 



(2/3 ~ 2/f ] 
1 - 2/3 



x+l _ L+l 

(2* - %f + 2 S ^ 



: lX )+ y^y? 

l<xi <X2<L 



1-2/3 

(2a - l)a 23 y? + (2a - l)yf 



_ L+l _ xi V X1+1 - li X2 

+ 2a^ ?L + 2aa 23 <7i 3 ^_J*L- + 2 flOJB .% % 



1-2/3 



1-2/3 



1-2/3 



|a;i,x 2 ) 



where 



Vj = eikpj . J = 1, 2, 3 and £j j£ = <r( yj , w ) , 1 < j ^ £ < 3. 



4 Let us stress the dependence on P € 63 in the definition of j/.,: it is used below. 
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After use of the Bethe equation, y% = 023CT13, it can be recasted as 



5+^3= E Ap M) E^ 1 ^ 3 ) 

Pee 3 I x=i 



2s . 

2S - 2 + (y 3 - (T230-13, 

1-2/3 



\x, x) 



+ E yiy? 

l<xi<X2<X/ 



<T2ZVT 2s - 1 + 2s 



2/3 



2.s 



+ yf 2s - 1 + 2s 



2s 



1-2/3 1-2/3 



0"23 



1-2/3 1-2/3 



£i,x 2 ) 



0"13 



Using the sum on P to relabel the variables yj, one can rewrite this equality as 



S+*3= £ ^p(k)|^i( M3 ) a 
PS6 3 I z=l 



(2s - 2)(l + CJ12 + (723 + C12C13 + C23<713 + <7120'13C'"23) 



2s 



+ ■: (2/3(1 — (T12) - <T23(Ti3 — <Tl2(Tl30'23) 

1-2/3 



2.s 



+ ~ (2/2<723(l - Cl3) - 0"12 - <7l2<7l3) 

1-2/2 



+ 



2.9 



1-2/1 



(yi 0-120-13(1 - o- 23 ) - 1 - 0-23) 



\x, x) 



+ E tar^ 2 

l<xi<3;2<L 



(1 + o-i 2 )(2s - 1) + 2s 



2/2 



1-2/2 1-2/1 



+ 2soi2 



2/1 



1 



l-2/i 1-2/2 



1^1, £2 



+ E yTiy^T 

l<xi<X2<L 



(1 + o- 23 )(2s - 1) + 2s 



2/3 



1-2/3 1-2/2 



+ 2S023 



2/2 



1-2/2 1-2/3 



|£l,Z2 



The term inside the square bracket [• • • ] in factor of (2/1 2/2 ) Xl 2/3 2 on the- one hand, and in factor 
of 2/1 1 (2/22/3) 312 011 the other hand, identically vanishes. It is in fact the same identity as the one 
used to show that 5 +x I / 2 = 0. It is also the identity used by Gaudin [12] to prove, for spin i, 
that Vm, is a highest weight vector. 

When the spin is higher than ^, it remains the term in factor of \x,x), which is a state that 
does not exist when s = |. The square bracket in front of \x, x) also identically vanishes, another 
identity due to the form of the scattering matrix a, and we get 5 + \I'3 = 0. 

When s = 1 this new identity is sufficient (together with the one used for S + ^2) to prove 
that ty m , Vm, is a highest weight vector. However, for s > 1, to prove that S + ^^ = 0, one needs 
to consider the state \x, x, x), that we will lead to another identity of the scattering matrix, and 
so on: for spin s, one needs 2s identities to prove that \& m , Vm, is a highest weight vector. 
Hence the difficulty to get a generic proof of it. 
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4 Conclusions 

In previous studies, the eigenfunctions of the spin s chain studied in this paper were known 
thanks to the algebraic Bethe ansatz. This later construction allows one to compute the cor- 
relation functions [171 118} \19\. Prior to that computation, the coordinate Bethe ansatz allowed 
Gaudin [20] to show, for spin k chains, orthogonality relations for the Bethe eigenfunctions, and 
prove a closure property for these functions. The explicit form of the eigenfunctions computed 
in this note is a first step toward a generalisation to spin s chains. The same method can also 
be applied to spin chains associated to higher rank algebras, for which less is known. 

In the same way, the spin 1 XXX chain with open (diagonal) boundaries has been studied 
in [21]: there is no doubt that their results can be generalized to spin s, using the present 
approach. The advantage of this method lies in the fact that we do not need to solve the reflection 
equation before computing the spectrum. We may start with general boundary conditions and 
find the ones for which the method is still consistent. In this way, the boundaries which keep 
the model solvable are classified. 

We also believe that the method presented here can be applied to solve the XXZ model with 
higher spin in the case of periodic boundary conditions. These cases have been treated through 
algebraic Bethe ansatz, see e.g. [22} 123] . More interestingly, general XXZ models with open 
boundary conditions can also be treated in this way, see e.g. [23] where a first account has been 
given. 

To conclude, we hope that this paper convinced the reader that the coordinate Bethe ansatz 
is a very powerful method and can be applied to solve a rather large class of integrable models. 
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